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(54) Precede de filtrage frequence! applique au debruitage de signaux sonores mettant en oeuvre 
un filtre de Wiener 



(57) Le procede selon I'invention utilise le filtrage 
frequentiel de Wiener pour debruiter des signaux sono- 
res (u(0). II comprend une etape preliminaire consistant 
a numeriser par echantillonnage des signaux sonores 
a debruiter (u(f)) et a les decouper en frames. Le pro- 
cede comprend ensuite une premiere serie d'etapes 
consistant a creer un modele de bruit (1) sur N trames, 
estimer (2) la densite spectrale du bruit, I'energie du mo- 



dele de bruit et calculer un coefficient traduisant la dis- 
persion statistique du bruit. II comprend egalement une 
seconde serie d'etapes consistant a calculer (3) la den- 
site spectrale des signaux a debruiter, pour chaque tra- 
me. Les coefficients du filtre de Wiener sont modifies 
(4), pour chaque trame successivement traitee, a I'aide 
des parametres determines a ('issue des deux series 
d'etapes, de maniere a introduire une compensation 
energetique et une surestimatbn adaptative du bruit. 
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Description 

[0001] La presente invention concerne un procede de filtrage frequentiel mettant en oeuvre un filtre de Wiener. 
[0002] Elle s'applique notamment, bien que non exclusivement, au debruitage de signaux sonores contenant de la 
s parole captee en milieux bruites et de facon plus generale au debruitage de tous signaux sonores. 

[0003] Les domaines principaux concernent les communications telephoniques ou radiotelephoniques, la reconnais- 
sance vocale, la prise de son a bord d'aeronefs civils ou militaires, et de facon plus generale de tous vehicules bruyants, 
les intercommunications de bord, etc. 

[0004] A titre d'exemple non limitatif, dans le cas d'un aeronef, les bruits resuitent des moteurs, de la climatisation, 
io de la ventilation des equipements de bord ou des bruits aerodynamiques. Tous ces bruits sont captes, au moins par- 
tiellement, par le microphone dans lequel parle !e pilote ou un autre membre de I'equipage. En outre, pour ce type 
d'application en particulier, une des caracteristiques des bruits est d'etre tres variables dans le temps. En effet, its sont 
tres dependants du regime de fonctionnement des moteurs (phase de decollage, regime stabilise, etc.). Les signaux 
utiles, c'est-a-dire les signaux representant les conversations, presentent egalement des particularites : ils sont le plus 
is souvent de breve duree. 

[0005] Entin, quelle que soit I'application envisaged, si on s'interesse au "voisement", on peut mettre en evidence 
certaines particularites. Comme il est ccnnu, le voisement concerne des caracteristiques elementaires de morceaux 
de parole, et plus precisement concerne les voyelles, ainsi qu'une partie des consonnes : "b", "d", "g", "j", etc. Ces 
lettres se caract6risent par un signal audiophonique de structure pseudo-p6riodique. 
20 [0006] En traitement de la parole, il est courant de considerer que les regimes stationnaires, notamment le voisement 
precite, s'etablissent sur des durees comprises entre 10 et 20 ms. Cet intervalle de temps est caracteristique des 
phenomenes elementaires de la production de la parole et sera denomme trame ci-apres. 

[0007] Aussi, il est usuel que les procede de debruitage prennent en compte cette caracteristique importante des 
signaux sonores comprenant de la parole. 

2S [0008] Ces precedes comprennent generalement les etapes principales suivantes : un decoupage en trames du 
signal audiophonique a debruiter, le traitement de ces trames par une operation de transformed de Fourier (ou d'une 
transformed similaire) pour passer dans le domaine frequentiel, le traitement de debruitage proprement dit par filtrage 
num6rique, et un traitement, dual du premier, par une transformed de Fourier inverse, pour revenir dans le domaine 
temporel. La derniere etape consiste en une reconstruction du signal. Cette reconstruction peut etre obtenue en mul- 

30 tipliant chacune des trames par une fenetre de ponderation. 

[0009] Un des filtres numeriques les plus utilises pour ce type d'application est le filtre de Wiener, en particulier un 
filtre de Wiener dit optimal. Celui-ci presente I'avantage de traiter de facon differenciee les trames successives. 
[001 0] En d'autres termes, et de facon plus generale, Le filtrage optimal de Wiener se trouve au centre des methodes 
optimales de traitement du signal, basees sur les caracteristiques statistiques du second ordre et done de la notion 

35 de correlation. 

[001 1 ] Le filtrage de Wiener permet la separation des signaux par decorrelation. Son importance est lied a la simplicity 
des calculs theoriques. En outre, il peut s'appliquer a une multitude de processus particuliers, et notamment, en ce 
qui concerne I'application pr6fered vised par I'invention, ('extraction d'un bruit polluant un signal de parole. 
[0012] Cependant, dans Tart connu, un probleme classique rencontre lors du debruitage par filtrage de Wiener est 

40 la presence d'un bruit, appele bruit musical, qui degrade la perception du signal debruite. Ce bruit musical est du aux 
fluctuations des densites spectrales du bruit present dans le signal d'entree. Pour certaines trames, en effet, la densite 
spectrale du bruit est superieure, au mcins sur un canal frequentiel, a celle du modele de bruit que Con utilise dans 
ces techniques. Dans ce cas, les mecanismes propres au filtrage de Wiener provoquent ('apparition d'un bruit residue! 
sur le signal debruite. Celui-ci est particulierement desagreable d'un point de vue perceptuel de part son instabilite. 

45 En effet, lors de Tecoute d'un signal de parole, on distingue des bruits residuels sous la forme de "gtouglou", qui 
s'apparente a des distorsions que Ton peut attribuer a une grande variabilite du bruit polluant le signal de parole debruite 
ou signal "utile". 

[0013] L'invention se fixe done pour but de pallier les inconvenients des precedes de filtrage de I'art connu : notam- 
ment I'inconvedient principal qui vient d'etre rappele : la presence d'un bruit residuel parasite dans le signal debruite, 
so dit "bruit musical". L'invention vise, de facon plus generale, a augmenter I'intelligibilite de la parole, dans son application 
principale. 

[0014] En vue d'attenuer fortement les effets du bruit musical, ('invention tire profit des deux observations experi- 
mentales suivantes : 

55 . la probability de bruit musical est d'autant plus forte que I'estimee des densites spectrales du bruit est instable 
d'une trame a I'autre ; 

la probability de presence de bruit musical est d'autant plus forte que I'estimee de la densite spectrale du bruit est 



2 
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faible par rapport a sa densite spectrale reelle. 

[0015] Selon une caracteristique principale de I'invention, le filtre de Wiener utilise pour le filtrage numerique est 
modifie de facon optimisee en y introduisant un terme de compensation energetique visant a surestimer te niveau de 
* 5 bruit. En outre, ce terme de compensation est adaptatif. 

[001 6] L'invention a done pour objet un procede de filtrage f requentiel pour le debruitage de signaux sonores bruites 
constitues de signaux sonores dits utiles melanges a des signaux de bruit, le procede comprenant au moins une etape 
de decoupage desdits signaux sonores en une serie de trames identiques d'une longueur determinee et une etape de 
filtrage frequentiel a ('aide d'un filtre de Wiener, caracterise en ce qu'il comprend, en outre, les etapes suivantes : 

10 

elaboration a partir desdits signaux bruites d'un modele de bruit sur un nombre N determine desdites trames, N 
etant compris entre des bomes minimale et maximale predeterminees ; 

application d'une transformee de Fourier auxdites N trames ; 

15 

estimation, pour chaque trame dudit modele, de la densite spectrale de cette trame ; 

estimation de la densite spectrale moyenne dudit modele de bruit ; 

20 - calcul, a partir de ces deux estimations, d'un coefficient de surestimation statistique, ledit coefficient statistique 
§tant egal au rapport maximal, pour lesdites A/trames du modele de bruit, entre le maximum de la densite spectrale 
d'une trame consideree dudit modele de bruit, et le maximum de la densite spectrale estimee du modele de bruit ; 

estimation, pour chaque trame desdits signaux a debruiter, de sa densite spectrale ; et 

25 

modification, pour chaque trame desdits signaux a debruiter, des coefficients dudit filtre de Wiener pour que la 
relation suivante soit verifiee : 



30 



1 - a • maxi 



Vv>J 



relation dans laquelle a et p sont des coefficients fixes predetermines, dits coefficient stattque de compen- 
sation energetique et coefficient d'attenuation exponentielle, respectivement, v decrit ('ensemble des canaux fre- 
35 quentiels de iadite transformee de Fourier, yJv) etant I'estimee de la densite spectrale de la trame a debruiter, y x 

(v) est Iadite densite spectrale du modele de bruit, et maxi ledit coefficient de surestimation statistique, modifiant 
le coefficient statique de compensation energetique a. 

[0017] L'invention sera mieux comprise et d'autres caracteristiques et avantages apparaltront a la lecture de la des- 
40 cription qui suit en reference aux figures annexees, parmi lesquelles : 



la figure 1 illustre, sous forme de bloc diagramme, les principales etapes du procede selon l'invention ; 
la figure 2 illustre schematiquement un filtre de Wiener de Tart connu ; 

45 

la figure 3 est un diagramme illustrant la densite spectrale d'un modele de bruit et les densites spectrales y u de 
chaque trame de ce modele de bruit ; 



les figures 4a et 4b sont des diagrammes comparatifs illustrant ces memes parametres avec surestimation de la 
so densite spectrale du modele de bruit ; 



la figure 5 est un diagramme illustrant ces memes parametres avec surestimation adaptative de la densite spectrale 
du modele de bruit ; 

ss . la figure 6 represente un exemple typique de signal issu d'une prise, de son bruitee ; 

la figure 7 est un organ igramme representant les etapes d*un procede particulier de recherche d'un modele de 
bruit ; 
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et la figure 8 est un organigramme detaille representant les etapes du procede de filtrage numerique selon un 
mode de realisation prefere de ('invention. 

[0018] Les principals phases et etapes du procede selon I'invention vont maintenant etre decrites par reference au 
s bloc diagramme de la figure 1 . Chaque bloc, references 0 a 5, represente une phase du procede, elle-meme pouvant 
etre subdivisee en etapes elementaires. 

[0019] Dans ce qui suit, pour fixer les idees et sans que cela limite en quoi que ce soit la portee de invention, on 
va se placer, dans le cadre du traitement de la parole bruitee. Comme il a ete indique precedemment, il est courant 
de considerer que les regimes stat ion na ires, notamment le voisement, s'etablissent sur des durees comprises entre 
io 10 et 20 ms, intervalle de temps caracteristique des phenomenes elementaires de la production de la parole et qui 
sera denomme trame ci-apres. 

[0020] Comme dans Tart connu, le precede de I'invention, comprend une etape de decoupage en trames du signal 
audiophonique a deb ru iter (bloc 0). 

[0021] Dans la pratique, on met en oeuvre des techniques numeriques. Aussi, les signaux de trame ne sont pas des 
'5 signaux a "evolution continue", mais des signaux discrets, obtenus par echantillonnage. On suppose que les signaux 
sont echantillonnes a la periode T e , avant traitement numerique. II est courant de considerer alors 2? echantillons pour 
une trame de signal, en choisissant p de maniere a ce que la valeur 2PT e soit de I'ordre grandeur de la duree D d'une 
trame. A titre d'exemple, pour une frequence d'echantillonnage de 10 kHz, on choisit souvent des trames de 12,8 ms, 
de maniere a pouvoir disposer de 128 points pour chaque trame, ce qui constitue une puissance de deux. Le nombre 
20 d'echantillons correspondant a une trame sera note ci-apres LGtrame. La relation suivante : D = LGtrame x T e est 
done satisfaite. L'etape de decoupage en trames, comme indique sur la figure 1 , est done pr^cedee d'une etape de 
numerisation par echantillonnage. 

[0022] Par convention, le signal d'entree sera note u(t) t le signal utile s(t) et le bruit perturbateur x(t) de telle facon 
que : 

25 

u(t) - s(t) + x(t) en temps continu (1 ) 

3Q u(/cT e ) = s( kT e ) + x{kT e ) en temps discret (2) 

[0023] Les etapes de numerisation et de decoupage en trames (bloc 0) sont communes a Tart connu. Les echantillons 
numeriques ainsi crees sont ranges dans une memoire tampon circulante de type "FIFO" (e'est-a-dire du type "premier 
entre - premier sorti") afin d'etre lus sous forme de trames successives. 
35 [0024] Les trames success ivement lues subissent alors une serie d'etapes de traitement autonomes, selon deux 
votes que Ton peut qualifier de a parallefes". 

[0025] Les operations effectuees dans le bloc 1 , consiste a identifier des segments du signal a debruiter ne contenant 
que du bruit. La sortie de ce bloc est constitute d'une suite d'echantillons numeriques representatifs du bruit seul. En 
d'autres termes, un modele de bruit est elabore a partir des signaux bruites, ou plus precisement a partir des trames 
40 successivement lues (bloc 0). De nombreux procedes peuvent etre mis en oeuvre et un exemple de procede de re- 
cherche de modele de bruit sera expiicite ci-apres. 

[0026] Dans le bloc 2, trois etapes sont realisees et consistent, a partir des echantillons fournis par le bloc 1, a 
effectuer : 

45 - I'estimation de la densite spectrale moyenne du bruit (par exemple par spectre moyen et correlogramme lisse) ; 
la determination de Penergie moyenne du modele de bruit ; 
et la determination d'un coefficient traduisant la dispersion statistique du bruit. 

so 

[0027] Les etapes ci-dessus, el notamment la derniere etape qui constitue une des caracteristiques principales de 
I'invention, seront detaillees ci-apres. 

[0028] Dans la branche "parallele", le bloc 3 comporie une etape d'estimation de la densite spectrale de la trame 
courante de signal et de calcul de son energie. 
55 [0029] Dans le bloc A t selon une autre caracteristique essentielle de i'invention, les coefficients du filtre frequentiel 
effectuant le debruitage du signal sont determines de la maniere qui sera detaillee ci-apres. Comme il a ete indique, 
le procede de I'invention est base sur une compensation 6nergetique et une surestimation du bruit. 
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[0030] Enfin, dans ie bloc 5, le signal temporel debruite est reconstruit, en assurant la meilleure continuity possible 
entre les trames. Dans d'autres applications que I'application principale visee par invention les signaux peuvent etre 
exploites tels quels par des divers precedes tels que la reconnaissance automattque de la parole. En soi, cette phase 
du procede est commune a Tart connu, et il n'y a pas lieu de detainer la methode de reconstruction ou d'exploitation 
des signaux en sortie du bloc 4. 

[0031] Selon la caracteristique principale de I'invention, le procede permet de modifier et d'optimiser les coefficients 
du filtre de Wiener utilise pour la phase de debruitage proprement dite (bloc 4), de facon a eliminer ou, pour le moins, 
fortement attenuer, les bruits parasites dits "musicaux". 

[0032] Comme il a ete rappeie, ces bruits sont attribuables a deux causes principales : 

a/ la probability de bruit musical est d'autant plus forte que I'estimee des densites spectrales du bruit est instable 
d'une trame a I'autre ; 

b/ la probability de presence de bruit musical est d'autant plus forte que I'estimee de la densite spectrale du bruit 
is est faible par rapport a la densite spectrale reelle du bruit 

[0033] Selon ['invention, en relation avec la cause a/, la dispersion est quantifiee par un coefficient issu de I'analyse 
effectuee dans le bloc 2, a partir du modele de bruit 6labore dans le bloc 1 . 

[0034] De m§me, en relation avec la cause b/, pour recluire ('influence de la densite spectrale du bruit, en particulier 
20 lorsqu'elle est faible, le procede selon I'invention effectue une surestimation de cette densite spectrale, en y introduisant 
un degre d'adaptivite afin d'optimiser la perception du signal debruite. 

[0035] Avant de decrrre plus en detail le procede de Pinvention, il est utile de rappeler brievement les caracteristiques 
d'un filtre de Wiener selon Tart connu. 

[0036] La figure 2 illustre de facon tres schematique un filtre de Wiener utilise pour debruiter un signal bruite U(n). 
25 [0037] A titre d'exemples non limitatifs, des filtres de Wiener sont decrits dans les livres suivants, auxquels on pourra 
se referer avec profit : 

Yves THOMAS : "Signaux et systemes lineaires', editions MASSON (1994) ; et : 

30 . Francois MICHAUT : "Methodes adaptatives pour le signal', Edition HERMES (1 992). 

[0038] Sur la figure 2 les conventions suivantes ont ete adoptees : 

U(n) : transformed de Fourier discrete du processus aleatoire observe, soit le signal bruite ; 

35 

S(n) : transformed de Fourier discrete du processus "desire - , a estimer par filtrage lineaire de U(n) ; 
X(n) : transformed de Fourier discrete du bruit additif polluant le signal utile ; 

A 

40 - SfnJ : estimation de S(n) exprimee dans le domaine de 

A 

Fourier, avec e=S-S= erreur d'estimation (S etant le signal debruite reel) ; et 
W(z) : filtre d'estimation exprime dans le domaine frequentiel. 



45 



SO 



55 



[0039] Le filtre optimal de Wiener minimise la distance entre les variables aleatoires S(n) et S(n) mesured par I'erreur 
quadratique moyenne J : 

J=E[(S(n)-Sfn;) 2 ) (3) 

[0040] La minimisation de ce critere revient a rendre I'erreur d'estimation orthogonale au signal observe, ce qui se 
traduit par le principe d'orthogonalrte : 

E[e(n).U*(n)] = 0 (4) 
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[0041] Ennotant: 

7 S la densite spectrale du signal utile, et 
Yx la densite spectrale du bruit parasite, 

le filtre de Wiener est decrit par la relation suivante : 

[0042] En prenant en compte I'independance de S^et de X(n), on obtient la relation ci-dessous : 

Yu=Ys + V x (6) ' 

relation dans laquelle y u representant la densite spectrale du signal observe. 
[0043] La relation decnvant le filtre de Wiener devient done finalement : 

y s ( n ) * yx( n ) m 

[0044] En pratique, e'est cette seconde formulation du filtre de Wiener qui est utilisee, puisqu'elle ne fait intervenir 
25 que des termes directement accessibles. c'est-a-dire : d'une part, le signal bruite recu du bloc 3 et, d'autre part, le bruit, 
prealablement determine par le calcul du modele de bruit (bloc 1). 

[0045] II doit etre remarque que les coefficients W(n) du filtre de Wiener sont toujours positifs. Si des artefacts de 
calcul provoque une valeur negative pour un coefficient, ce coefficient est rendu egal a zero. 

[0046] Selon Tart connu, la suppression du bruit additif par une methode de soustraction spectrale, telle qu'elle est 
30 realisee par un filtre Wiener, debouche sur la creation de bruits dits "musicauxV Pour eviter I'apparition de ces bruits 
parasites desagreables a I'ecoute et nuisibles a rintelligibilite de la parole, ou pour le moins empecher au maximum 
leur apparition, selon une caracteristique essentielle de I'invention, les coefficients du filtre de Wiener sont modifies a 
I'aide de parametres determines dans les blocs 2 et 3, de la maniere qui va maintenant etre detaillee. 
[0047] Lorsque le signal cfentree ne contient que du bruit, le 'bruit musical' supplemental est present parce que, 
dans la pratique, I'estimation du rapport ^ fluctue a chaque frequence, bien qu'en theorie ce rapport devrait etre egal 
a I'unite quelles que soient les frequence^. Ce sont ces erreurs d'estimation qui produisent des filtres attenuateurs 
dont les variations des coefficients sont aleatoires, selon les frequences et au cours du temps. 
[0048] Pour fixer les idees, on considere I'exemple du debruitage d'un bruit seul, echantillonne a 44 kHz. On deter- 
mine la densite spectrale y x d'un modele de bruit choisi a Paide de ce signal et les densites spectrales y a de chaque 
40 trame (de longueur LGtrame) de ce bruit. 

[0049] On a represents la variation de ces deux parametres sous forme de courbes dans le dagramme de la figure 
3, en fonction du nombre de canaux de transfomnee de Fourier FFT. Pour tracer les courbes, il a ete suppose que la 
longueur de trame etait de 128 echantillons, soit LGlrame=123. 

[0050] Ce diagramme montre clairement que les allures des deux courbes y x et y 0 sont similaires mais les deux 
45 estimees presented une difference cfamplitude nette. Le pic principal de y v qui se situe a la frequence 2.75 kHz (64 
canaux FFT correspondant a 22 kHz, soit la demi-frequence d'echantillonnage) a une amplitude environ sept fois 
superieure a celui de y x situe a la meme frequence. Ceci constitue la raison principale de la presence des bruits 
"musicaux": Lorsque, pour certaines frequences ret erencees v, yjv) est bien superieur a y^v), cela signifie, en theone, 
que la trame ne contient pas seulement du bruit mais une autre partie de signal. Dans ce cas, le filtrage de Wiener 
so selon I'art connu debruite la trame correspondante comme si elle contenait du signal de parole utile, ce qui entraine 
la presence de residus de bruits. 

[0051] Pour eviter cet effet parasite, le precede selon I'invention modifie de facon optimisee les coefficients du filtre 
de Wiener et intrbduit un terme de compensation energetique, venant surestimer artificiellement le niveau du bruit, 
avec differents niveaux cfadaptativite de cette compensation. 
55 [0052] Les coefficients du filtre de Wiener mocfme obeissent a la relation suivante : 
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W(v) = [ 1 - a • maxi • ~ - ^4 | ( 8 ) 

- 5 [0053] En se reportant de nouveau a la relation (7), on constate aisement que quatre nouveaux termes ont ete 
introduits, a savoir : 

p : coefficient d'attenuation exponentielle ; 

10 a : coefficient statique de compensation energetique ; 

rapport de ponderation energetique ; et 

u 

maxi. coefficient de surestimation statistique issu de I'analyse statistique du bruit, ce a partir d'un modele de bruit 
is etabii lors de la phase du procede correspondant au bioc 1. 

[0054] Chacun de ces termes va maintenant etre explicite. 

[0055] Le coefficient d'attenuation exponentielle p est un terme communement utilise dans la litterature consacre 
au dornaine du filtrage numerique et, plus particulierement, au d6bruitage. Une valeur typique de ce parametre est de 
20 o,5. 

[0056] A titre d'exemple non limitatif, on pourra se reporter a Particle de L. Arslan : A. Mc Cree et V Viswana- Than, 
intitule :"New Methods for adaptive noise suppression", IEEEE, mai 1995, pages 812-815. 

[0057] Le coefficient de compensation energetique statique a permet de surestimer le bruit et est particulierement 
pertinent dans le cas de la suppression de bruit seul. En effet, une valeur typique de <x=10 appliquee a Pexemple de 
25 la figure 3 augmente I'estimeedu spectre moyen de bruit y x d'environ +10 dB, ce qui permet alorsdediminuer le niveau 
de bruit residuel, puisque les coefficients du filtre de Wiener ne peuvent etre n6gatifs. Dans le cas contraire, Ms sont 
alors forces a zero. 

[0058] Cependant, si cette modification est tres efficace pour eliminer le bruit seul, elle pose a son tour des problemes 
lorsque les trames a debruiter contiennent du signal utile. Si ce signal utile est beaucoup plus energetique que le bruit, 
30 ce coefficient multipiicateur a n'a pas d' effet sur la degradation de ce signal. Mais, dans le cas contraire, il peut exister 
des frequences v pour lesquelles une trame de signal utile a une energie non negligeable mais proche de celle du 
bruit pour les memes frequences. Dans ce cas, la multiplication par a de y^v) impose des coefficients de Wiener W 
(v) nuls et done entraTne une disparition de I'energie du signal pour ces frequences. 

[0059] Ce probleme est illustre par les figures 4a et 4b. Sur ces figures les conventions suivantes ont ete adoptees. 

35 

y u : : densite spectrale de la trame de signal consideree (trame de signal faiblement energetique devant le bruit) ; et 
y x : densite spectrale du modele de bruit choisi (bloc 1 ). 

[0060] La courbe de la figure 4a permet de constater que I'energie du signal dans la bande de frequences Av, re- 
40 presentee par la densite spectrale y^ n'est pas negligeable. 

[0061] En se referant a la figure 4b, on peut constater que la multiplication de y x par le parametre ct=10 rend a.y* 
superieura Y„dans la bande Av. II s'ensuit que le gain de Wiener est nul pour cette bande de frequences qui n'apparalt 
plus dans la trame debruitee. 

[0062] Le rapport de ponderation energetique decrit ci-dessous permet de reduire cette distorsion dans le signal 
45 d^bruite. 

[0063] Comme indique precedemment, le debrurtage du bruit seul est correct, mais il peut etre trop brutal dans les 
parties du signal utile. 

[0064] Dans une variante preleree de Tinvention, on remedie a cet inconvenient en faisant varier le coefficient a, ce 
en fonction de la presence ou non d'une partie de signal utile dans le signal a debruiter. De facon avantageuse, a reste 

so proche de d'une valeur typique egale a 10, lorsque le signal bruite ne contient que du bruit, et varie entre 0 et 10, 
lorsqu'un signal utile est present dans le signal bruite. E On introduit done avantageusement un degre d'adaptativite. 
[0065] C'est la fonction qui est assignee au rapport qui vient multiplier a dans la relation (8), rapport dans lequel 
£ x est I'energie moyenne du modele de brurt et E w I'enetgie de la trame courante. Cela permet done aux coefficients 
du filtre de Wiener de changer a chaque trame de facon differenciee selon la presence plus ou moins grande (en terme 

ss d'energie) du signal de parole. 

[0066] Si E x == Ey, alors a=10 et la trame est consideree comme du bruit seul. Elle est correctement debruitee. 
[0067] Si au contraire E x « E^ cela signrfie que la trame consideree est tres energetique devant le bruit et qu'il est 
necessaire d'attenuer au minimum cette partie de signal. 
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[0068] Cette troisieme modification est tllustree par la figure 5. Sur cette figure, la trame de signal consideree est la 
meme que celle utilisee pour les figures 4a et 4b, 

E x 

s a = 10et=^ = 0,2. 

fc u 

[0069] Grace a cette ponderation du coefficient a par-^sc, la bande de frequences Av' dans laquelle le signal utile 

est elimine, (c'est-a-dire les frequences pour lesquelles T« coefficients de y x sont superieurs a ceux de yj est bien 
10 moins importante que lors de la modification par multiplication du seul coefficient a=10. 

[0070] Ce type de filtre presente done une bonne efficacite en termes d'elimination des segments de signal degrades 

dans lesquels la parole est absente et de diminution des distorsions infligees au signal de parole utile. 

[0071] La probability de generation du "bruit musical" est egalement liee, comme il a ete indique, a la variance des 

estimees de la densite spectrale du bruit sur I'ensemble des trames. 
is [0072] En effet, plus les densites spectrales estimees du bruit varient d'une trame a I'autre, plus la formation du bruit 

"musical' 1 est probable. 

[0073] Selon un autre aspect important de I'invention, on rend dependant la valeur du coefficient de surestimation 
des proprietes statistiques du bruit. Pour ce faire, il est introduit un coefficient, appele max/'ci-apres, proportionnel a 
la dispersion des valeurs de densites spectrales du bruit. 
20 [0074] Le coefficient de surestimation devient alors : 

a = a* max/, avec maxi satisfaisant la relation suivante : 



25 

maxi = Max 

30 



relation dans laquelle : 

35 - N est le nombre de trames du modele de bruit ; 

v decrit I'ensemble des canaux frequentiels, soit LGtrame/2 canaux ; 

yfrjesi la densite spectrale de la i erne trame du modele de bruit dans le canal v ; et 

y x (v) est la densite spectrale du modele de bruit. 

40 [0075] Le coefficient max/ est egal au rapport maximal, pour toutes les trames du modele de bruit, entre le maximum 
de la densite spectrale de la trame du modele de bruit consideree, et le maximum de la densite spectrale estimee du 
modele de bruit. 

[0076] En d'autres termes, ce coefficient caractertse la disparite maximale du bruit pour les canaux frequentiels 
portant une energie importante. Multiplie par le coefficient ct : il apporte une attenuation complementaire proportionnelle 
45 a cette disparite. 

[0077] Pour elaborer une partie des parametres entrant dans la modification des coefficients du filtre de Wiener, il 
est necessaire de disposer d'un modele de bruit (bloc 1 de la figure 1). 

[0078] L'elaboration d'un modele de bruit d'un signal bruite est une operation classique en soi. Cependant, la methode 
specifique mise en oeuvre pour cette operation peut etre une methode de Tart connu, mais aussi une methode originale. 
so [0079] On va decrire ci-apres, par reference aux figures 6 et 7, une methode d'elaboration d'un modele de bruit, 
particulierement adaptee aux applications principales visees par le precede de I'invention, notamment le debruitage 
de signaux de parole bruites. 

[0080] La methode repose sur une recherche permanente et automatique d'un modele de bruit. Cette recherche est 
faite sur les echantitlons de signal u(t) numerises et stockes dans une memoire tampon d'entree. Cette memoire est 
55 capable de memoriser simultanement tous les echantillons de plusieurs trames du signal d'entree (au moins 2 trames 
et, dans le cas general, N trames). 

[0081] Le modele de bruit recherche est constitue par une succession de plusieurs trames dont la stabilite en energie 
et le niveau d'energie relative font penser qu'il s'agit d'un bruit ambiant et non d'un signal de parole ou d'un autre bruit 
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perturbateur. On verra plus loin comment se fait cette recherche automatique. 

[0082] Lorsqu'un modele de bruit est trouve, tous les echantillons des A/trames successives representant ce modele 
de bruit sont conserves en memoire, de sorte que le spectre de ce bruit peut etre analyse et peut servir au debruitage. 
Mais la recherche automatique de bruit continue a partir du signal d'entree u(t) pour trouver eventuellement un modele 
s plus recent et plus adapte, sort parce qu'il represente mieux le bruit ambiant, soit parce que le bruit ambiant a evolue. 
Le modele de bruit plus recent est mis en memoire a la place du precedent, si la comparaison avec le precedent montre 
qu'il est plus representatif du bruit ambiant. 

[0083] Les postulats de depart pour I'elaboration automatique d'un modele de bruit sont les suivants : 

io - le bruit qu'on veut eliminer est le bruit de fond ambiant, 

le bruit ambiant a une energie relativement stable a court terme, 

la parole est le plus souvent precedee d'un bruit de respiration du pilote qu'il ne faut pas confondre avec le bruit 
ambiant; mais ce bruit de respiration s'eteint quelques centaines de millisecondes avant la premiere emission de 
parole proprement drte, de sorte qu'on ne retrouve que le bruit ambiant juste avant ('emission de parole, 
is - el cntm. les bruits et la parole se superposent en termes d'energie de signal, de sorte qu'un signal contenant de 
in oarole ou un bruit perturbateur, y compris la respiration dans le microphone, contient forcement plus d'energie 
qu'un signal de bruit ambiant - ' 

— \OOS4\- ikjcn_/esulte qu'on fera I'hypothese simple suivante : le brurt ambiant est un signal presentant une energie 
20 mminviio sirtDic a court terme. Par court terme, il faut entendre quelques trames, et on verra dans I'exemple pratique 
dunnc ci-npieb que le nombre de trames destine a evaluer la stabilite du bruit est de 5 a 20. L'energie doit etre stable 
sur ptusicurs trames. laute de quoi on doit suppose r que le signal contient plutot de la parole ou un bruit autre que le 
bruit rtmbiant Elle doit etre minimale, faute de quoi on considere que le signal contient de la respiration ou des elements 
phonehques de parole ressemblant a du bruit mais se superposant au bruit ambiant. 
2$ [0085] La figure 6 represente une configuration typique devolution temporelle de l'energie d'un signal microphonique 
au moment d'un debut d'emission, de parole, avec une phase de bruit de respiration, qui s'eteint pendant quelques 
dizaines a centaines de millisecondes pour faire place au bruit ambiant seul, apres quoi un niveau d'energie eleve 
indique la presence de parole, pour revenir enfin au bruit ambiant. 

[0086] La recherche automatique du bruit ambiant consiste alors a trouver au moins N1 trames successives (par 
so exemple N1 = 5) dont les energies sont proches les unes des autres, e'est-a-dire que le rapport entre l'energie de 
signal contenue dans une trame et l'energie de signal contenue dans la ou, de preference, les trames precedentes 
est situe a Pinterieur d'une gamme de valeurs determinee (par exemple compris entre 1/3 et 3). Lorsqu'une telle suc- 
cession de trames d'energie relativement stable a ete trouvee, on stocke les valeurs numeriques de tous les echan- 
tillons de ces A/trames. Cet ensemble de NxP echantillons constitue le modele courant de bruit. II est utilise dans le 
35 debruitage. L'analyse des trames suWantes continue. Si on trouve une autre succession d'au moins N1 trames suc- 
cessives repondant aux memes conditions de stabilite d'energie (rapports d'energies de trames dans une gamme 
determinee), on compare alors l'energie moyenne de cette nouvelle succession de trames a l'energie moyenne du 
modele stocke, et on remplace ce dernier par la nouvelle succession si le rapport entre l'energie moyenne de la nouvelle 
succession et l'energie moyenne du modele stocke est inferieur a un seuil de remplacement determine qui peut etre 
40 de 1 ,5 par exemple. 

[0087] De ce remplacement d'un modele de bruit par un modele plus recent moins energetique ou pas beaucoup 
plus energetique, il resulte que le modele de bruit se cale globalement sur le bruit ambiant permanent. Meme avant 
une prise de parole, precedee d'une respiration, il existe une phase ou le bruit ambiant seul est present pendant une 
duree suffisante pour pouvoir etre pris en compte comme modele de bruit actif. Cette phase de bruit ambiant seul, 
45 apres respiration, est breve. Le nombre N1 est choisi relativement faible, afin qu'on ait le temps de recaler le modele 
de bruit sur le bruit ambiant apres la phase de respiration. 

[0088] Si le bruit ambiant evolue lentement, revolution sera prise en compte du fait que le seuil de comparaison 
avec le modele stocke est superieur a 1. S'il evolue plus rapidement dans le sens croissant, revolution risque de ne 
pas etre prise en compte, de sorte qu'il est preferable de prevoir de temps en temps une reinitialisation de la recherche 
so d'un modele de bruit. Par exemple, dans un avion au sol a I'arret, le bruit ambiant sera relativement faible, et il ne 
faudrait pas qu'au cours de la phase de decollage le modele de bruit reste fige sur ce qu'il etait a I'arret du fait qu'un 
modele de bruit n'est remplace que par un modele moins energetique ou pas beaucoup plus energetique. On expliquera 
plus loin les methodes de reinitialisation envisagees. 

[0089] La figure 7 represente un organigramme des operations de recherche automatique d'un modele de bruit 
55 ambiant. 

[0090] Le signal d'entree u (t) , echantiilonne a la frequence F e = 1/T e et numerise par un convertisseur analogique- 
numerique, est stocke dans une memoire tampon capable de stocker tous les echantillons d'au moins 2 trames. 
[0091] Le numero de la trame courante dans une operation de recherche d'un modele de bruit est designe par n et 
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est compte par un compteur au fur et a mesure de la recherche. A I'initialisation de la recherche, n est mis a 1. Ce 
numero nsera increments au fur et a mesure de ('elaboration d'un modele de plusieurs trames successives. Lorsqu'on 
analyse latramecourante n, le modele comprend deja par hypothese r>1 trames successives repondant aux conditions 
imposees pour faire partie d'un modele. 
5 [0092] On considere d'abord qu'il s'agit d'une premiere elaboration de modele, aucun autre modele precedent n'ayant 
ete construit. On verra ensuite ce qui se passe pour des elaborations ulterieures. 

[0093] L'energie de signal de la trame est calculee par sommation des carres des valeurs numeriques des echan- 
tillons de la trame. Elle est conservee en memoire. 

[0094] On lit ensuite la trame suivante de rang n = 2, et son energie est calculee de la meme maniere. Elle est 

io egalement conservee en memoire. 

[0095] On calcule le rapport entre les energies des deux trames. Si ce rapport est compris entre deux seuils S et S' 
dont I'un est superieur a 1 et I'autre est inferieur a 1 , on considere que les energies des deux trames sont proches et 
que les deux trames peuvent faire partie d'un modele de bruit. Les seuils S et S'sont de preference inverses I'un de 
I'autre (S = 1/S) de sorte qu'il suffit de definir I'un pour avoir I'autre. Par exemple, une vateur typique est S = 3, S = 

15 1/3. Si les trames peuvent faire partie d'un meme modele de bruit, les echantillons qui les composent sont stockes 
pour commencer a construire le modele, et la recherche continue par iteration en incremental n d'une unite. 
[0096] Si le rapport entre les energies des deux premieres trames sort de I'intervalle impose, les trames sont decla- 
res incompatibies et la recherche est reinitialisee en remettant na 1. 

[0097] Dans le cas ou la recherche continue, on incr6mente le rang n de la trame courante, et on effectue, dans une 
20 boucle de procedure iterative, un calcul d'energie de la trame suivante et une comparaison avec l'energie de la trame 
precedente ou des trames precedentes, en utilisant les seuils Set S\ 

[0098] On notera a ce propos que deux types de comparaison sont possibles pour ajouter une trame a n-1 trames 
precedentes qui ont deja ete considerees comme homogenes en energie : le premier type de comparaison consiste 
a comparer uniquement l'energie de la trame n a l'energie de la trame Le deuxieme type consiste a comparer 
25 l'energie de la trame n a chacune des trames 1 a n-1 . La deuxieme maniere aboutit a une plus grande homogeneite 
du modele mais elle a I'inconvenient de ne pas prendre en compte suffisamment bien les cas ou le niveau de bruit 
crott ou decroit rapidement. 

[0099] Ainsi, l'energie de la trame de rang n est comparee avec l'energie de la trame de rang n-1 et eventuellement 
d'autres trames precedentes (pas forcement toutes d'ailleurs). 
so [0100] Si la comparaison indique qu'il n'y a pas homogeneite avec les trames precedentes, du fait que le rapport 
des energies n'est pas compris entre 1/S et S, deux cas sont possibles : 

ou bien n est inferieur ou egal a un nombre minima! N1 en dessous duquel le modele ne peut pas etre considere 
comme significatif du bruit ambiant parce que la duree d'homog^neite est trop courte; par exemple N1 = 5; dans 
35 ce cas on abandonne le module en cours d'elaboration, et on reinitialise la recherche au debut en remettant n a 1 ; 

ou bien nest superieur au nombre minimal N1. Dans ce cas, puisqu'on trouve maintenant un manque d'homoge- 
neite, on considere qu'il y a peut-etre un debut de parole apres une phase de bruit homogene, et on conserve a 
titre de modele de bruit tous les echantillons des n-1 trames de bruit homog&nes qui ont precede le manque 
40 d'homogeneite. Ce modele reste stocke jusqu'a ce qu'on trouve un modele plus recent qui semble egalement 

representer du bruit ambiant. La recherche est reinitialisee de toute fa5on en remettant n a 1 . 

[0101] Mais la comparaison de la trame n avec les precedentes aurait pu encore aboutir a la constatation d'une 
trame encore homogene en 6nergie avec la ou les precedentes. Dans ce cas, ou bien n est inferieur a un deuxieme 
45 nombre N2 (par exemple N2- 20) qui represente la longueur maximale souhaitee pour le modele de bruit, ou bien n 
est devenu egal a ce nombre N2. Le nombre N2 est choisi de maniere a limiter le temps de calcul dans les operations 
ulterieures d'estimation de density spectrale de bruit. 

[0102] Si nest inferieur a N2, la trame homogene est ajouteeaux precedentes pour contribuer a construire le modele 
de bruit, n est increments et la trame suivante est analysee. 
so [0103] Si n est 6gal a A/2, la trame est egalement ajoutee aux n-1 trames homogenes precedentes et le modele de 
n trames homogenes est stocke pour servir dans Telimination du bruit. La recherche d'un modele est par ailleurs 
reinitialisee en remettant na 1. 

[0104] Les etapes precedentes concernent la premiere recherche de modele. Mais une fois qu'un modele a et6 
. stocke. il peut a tout moment etre remplace par un modele plus recent. 
55 [0105] La condition de remplacement est enccre une condition d'energie, mais cette fois elle porte sur l'energie 
moyenne du modele et non plus sur l'energie de chaque trame. 

[0106] Par consequent, si un modele possible vient d'etre trouve, avec A/trames ou N1 < N.< N2, on calcule l'energie 
moyenne de ce modele qui est la somme des energies des N frames, divisee par A/, et on la compare a l'energie 
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moyenne des A/ trames du modele precedemment stocke. 

[0107] Si le rapport entre I'energie moyenne du nouveau modele possible et I'energie moyenne du modele actuel 
en vigueur est inferieur a un seuil de remplacement SR, le nouveau modele est considere comme meilleur et on le 
stocke a la place du precedent. Sinon, le nouveau modele est rejete et I'ancien reste en vigueur. 

.5 [01 08] Le seuil SR est de preference legerement superieur a 1 . 

[0109] Si le seuil SR etait inferieur ou egal a 1 , on stockerait a chaque fois les trames homogenes les moins ener- 
getiques, ce qui correspond bien au fait qu'on considere que le bruit ambiant est le niveau d'energie au dessous duquel 

* on ne descend jamais. Mais, on eliminerait toute possibility devolution du modele si le bruit ambiant se mettait a 
augmenter. 

io [0110] Si le seuil SR etait trop eleve au dessus de 1 , on risquerait de mal distinguer le bruit ambiant et d'autres bruits 
perturbateurs (respiration), voire meme certains phonemes qui ressemblent a du bruit (consonnes sifflantes ou chuin- 
tantes par exemple). L'elimination de bruit a partir d'un modele de bruit cale sur la respiration ou sur des consonnes 
sifflantes ou chuintantes risquerait alors de nuire a I'intelligibilite du signal debruite. 

[0111] Dans un exemple prefere le seuil SR est d'environ 1 ,5. Au-dessus de ce seuil on conservera I'ancien modele ; 
is en dessous de ce seuil on remplacera I'ancien modele par le nouveau. Dans les deux cas, on reinitialisera la recherche 
en recommencant la lecture d'une premiere trame du signal d'entree u(t), et en mettant na 1 . 

[0112] Pour rendre ('elaboration du modele de bruit plus fiable, on peut prevoir que la recherche d'un modele est 
inhibee si une emission de parole est detectee dans le signal utile. Les traitements numeriques de signal couramment 
utilises en detection de parole permettent d'identifier la presence de paroles en se fondant sur les spectres caracte- 
20 ristiques de periodicite de certains phonemes, nolamment les phonemes correspondant a des voyelles ou a des con- 
sonnes voisees. 

[011 3] Le but de cette inhibition est d'eviter que certains sons soient pris pour du bruit, alors que ce sont des pho- 
nemes utiles, qu'un modele de bruit fonde surces sons soit stocke et que la suppression du bruit posterieure a I'ela- 
boration du modele tende alors a supprimer tous les sons similaires. 
25 [0114] Par ailleurs, il est souhaitable de prevoir de temps en temps une reinitialisation de la recherche du modele 
pour permettre une remise a jour du modele alors que les augmentations du bruit dmbiant n'ont pas ete prises en 
compte du fait que SR n'est pas beaucoup superieur a 1 . 

[0115] Le bruit ambiant peut en effet augmenter de facon importante et rapide, par exemple pendant la phase de- 
celeration des moteurs d'un avion ou d'un autre vehicule, aerien, terrestre ou maritime. Mais le seuil SR impose que 

so le modele de bruit precedent soit conserve lorsque I'energie moyenne de bruit augmente trop vite. 

[0116] Si on souhaite remedier a cette situation, on peut proceder de diffe>entes manieres, mais la maniere la plus 
simple est de reinitialiser le modele periodiquement en recherchant un nouveau modele et en I'imposant comme modele 
actif independamment de la comparaison entre ce modele et le modele precedemment stocke. La periodicite peut etre 
basee sur la duree moyenne d'elocution dans ('application envisagee ; par exemple les durees d'elocution sont en 

35 moyenne de quelques secondes pour Pequipage d'un avion, et la reinitialisation peut avoir lieu avec une periodicite 
de quelques secondes. 

[0117] La mise en oeuvre de la methode d'elaboration d'un modele de bruit (figure 1 : bloc 1) et, de facon plus 
generale du precede selon ['invention, peut se faire a partir de calculateurs non specialises, pourvus de programmes 
de calcui necessaires et recevant les echantillons de signaux numerises tels qu'ils sont fournis par un convertisseur 
^0 analogique-numerique, via un port adapte. 

[011 8] Cette mise en oeuvre peut aussi se faire a partir d'un calculateur specialise a base de processeurs de signaux 
numeriques, ce qui permet de traiter plus rapidement un plus grand nombre de signaux numeriques. 
[0119] Les calculateurs sont associes, comme il est bien connu, a differents types de memoires, statiques et dyna- 
mique, pour enregistrer les programmes et les doonees intermediates, ainsi qu'a des memoires circulantes de type 
•FIFO". Le systeme comprend enfin un convertisseur analogique-numerique, pour la numerisation des signaux u(t), 
et un convertisseur numerique-analogique, en tant que de besoin, si les signaux debruites doivent etre utilises sous 
forme analogique. 

[0120] En conclusion, et pour decrire de facon plus detaijlee le proedde de ('invention, on peut decouper les etapes 
de facon differente de ce qui a ete decrit en reference a la figure 1 (qui illustre le precede de facon plus synthetique). 
so La figure 8 est un diagramme resumant toutes les etapes du precede de filtrage selon ('invention, dans un mode de 
realisation prefere. 

[0121] Ces etapes se repartissent en un premier sous-ensemble d'etapes permettant de determiner les parametres 
dependant du modele de bruit, et un second sous-ensemble d'etapes permettant de determiner les parametre depen- 
dant seulement de la trame courante du signal a debruiter. 
55 [0122] La premiere etape du premier sous-ensemble, comprend une 6tape initiate de selection d'un modele de bruit 
adapte a ('application specifique, avantageusement un modele de bruit determine par la methode decrite ci-dessus, 
en reference aux figures 6 et 7. 

[0123] Ce premier sous-ensemble d'etapes comprend deux branches. 
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( [0124] Dans la premiere branche, on calcule pour chaque trame du modele de bruit (dans le domaine temporel). 

I'energie de la trame, puis on calcule I'energie moyenne des trames du modele, ce qui permet d'estimer I'energie 
moyenne du modele, c'est-a-dire le parametre E x 

[0125] Dans la deuxieme branche. on applique une transformee de Fourier aux trames du modele de bruit, de facon 
5 a passer dans le domaine frequentiel. Puis on determine successivement la densite spectrale de la trame / (avec / = 
1 . . N) du modele de bruit dans le canal frequentiel v, sort y/vj, et la densite spectrale du modele de bruit dans le canal 
frequentiel v, soit y x (v). A partir de ces deux parametres, on determine le coefficient statistique maxide telle sorte qu'il 
verifie la relation (9). Le parametre y x (v) est egalement utilise pour le calcul d'un des autres coefficients du filtre de 
Wiener. 

ro [0126] Le second sous-ensemble d'etapes comprend egalement deux branches. 

[0127] Dans la premiere branche, on determine I'energie de la trame courante, soit E^, et dans la seconde branche, 
on effectue I'estimation de la densite spectrale de la trame courante y w 

[0128] A partir de ces deux parametres et des parametres y x e\ E x , determines precedemment, on obtient les coef- 
licicnls [E./E J et [yjvj/yjv)). 

is [01 29] Tous les coefficients du filtre de Wiener, conforme a la relation (8), sont done determines a I'issu de ces etapes. 

Lcs coeticients a et p sont des coefficients fixes predetermines, typiquement egaux a 10 et 0,5, respectivement. 

[0130) A la lecture de ce qui precede, on constate aisement que ('invention atteint bien les buts qu'elle s'est fixes. 

[01 31 1 II doit etre clair cependant que invention n'est pas limitee aux seuls exemples de realisations explicitement 

deems notammcnt en relation avec les figures 1 a 8. 
20 [0132] En particulter, les exemples numeriques n'ont ete donnes que pour mieux preciser ('invention mais sont es- 

bcnticl emon: hc-b a I'application specifique envisagee. De ce fait, its participent d'un simple choix technologique a la 

porlco dc rHorrme du Metier. 

[01 33] En outre, comme il ete rappele, ('invention ne se reduit pas au seul domaine du filtrage de signaux contenant 
de la pcircle bruitee. meme si ce domaine constitue une des applications preferees. 

25 

Revendications 

1 . Precede de filtrage frequentiel pour le debruitage de signaux sonores bruites (u(t)) constitues de signaux sonores 
30 dits utiles melanges a des signaux de bruit, le procede comprenant au moins une etape de decoupage (0) desdits 

signaux sonores en une serie de trames identiques d'une longueur determinee et une etape de filtrage frequentiel 
(4) a I'aide d'un filtre de Wiener, caracterise en ce qu'il comprend, en outre, les etapes suivantes : 

elaboration a partir desdits signaux bruites (u(t)) d'un modele de bruit ( 1 ) sur un nombre N determine desdites 
55 trames, A/etant compris entre des bornes minimale et maximale predeterminees ; 

application d'une transformee de Fourier auxdites N trames ; 

estimation (2), pour chaque trame dudit modele, de la densite spectrale de cette trame ; 
estimation (2) de la densite spectrale moyenne dudit modele de bruit ; 
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calcul (2), a partir de ces deux estimations, d'un coefficient de surestimation statistique, ledit coefficient sta- 
tistique etant egal au rapport maximal, pour lesdites N trames du modele de bruit, entre le maximum de la 
45 densite spectrale d'une trame consideree dudit modele de bruit, et le maximum de la densite spectrale estimee 

du modele de bruit ; 

estimation (3), pour chaque trame desdits signaux a debruiter (u(t)), de sa densite spectrale ; et : 

so - modification (4), pour chaque trame desdits signaux a debruiter (u(t)) t des coefficients dudit filtre de Wiener 

pour que la relation sutvante soit verifiee : 
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LV(v) = | 1 -a.max/.^yj , 

relation dans laquelle a et P sont des coefficients fixes predetermines, dits coefficient statique de compensation 
energetique et coefficient d'attenuation exponentielle, respectivement, v decrit I'ensemble des canaux fre- 
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quentiels de ladite transformee de Fourier, y^vjetant I'estimee de la densite spectrale de ta trame a debruiter, 
yjv) est ladite densite spectrale du modele de bruit, et maxi ledit coefficient de surestimation statistique, 
modifiant le coefficient statique de compensation energetique a. 

Precede selon la revendication 1 , caracterise en ce que ledit coefficient statistique maxi verifie la relation suivante : 



Procede selon Tune des revendications 1 ou 2, caracterise en ce qu'il comprend les etapes supplementaires 
suivantes : 

calcul de I'energie moyenne dudit modele de bruit E x ; 

calcul, pour chaque trame desdits signaux a debruiter (u(t)) t de I'energie de la trame en cours E u ; et 

multiplication dudit coefficient statique de compensation energetique a par un coefficient de ponderation ener- 
getique egal au rapport EJE U de maniere a modifier selectivement ces coefficients pour chaque trame desdits 
signaux a debruiter (u(t)) par application d'un coefficient continument variable entre un extrema et un minima, 
I'extrema etant sensiblement egal a Tunite lorsque lesdits signaux utiles sont absents desdits signaux a de- 
bruiter (u(t)) et sensiblement egal a zero lorsque I'energie desdits signaux utiles est tres superieure a I'energie 
desdits signaux de bruit, et a ce que lesdits coefficients du filtre de Wiener satisfassent la relation suivante : 



Procede selon Tune quelconque des revendications precedentes, caracterise en ce que ledit coefficient statique 
de compensation energetique a est egal a 10. 

Procede selon Tune quelconque des revendications precedentes, caracterise en ce que ledit coefficient d'attenua- 
tion exponentielle p est egal a 0,5. 

Procede selon Tune quelconque des revendications precedentes, caracterise en ce qu'il comprend une etape 
initiale (0) consistant a numeriser lesdits signaux a debruiter (u(t)) par echantillonnage, chaque trame comprenant 
p echantillons. 

Procede selon la revendication 6, caracterise en ce que ledit modele de bruit (1) est obtenu par une recherche 
repetitive effectuee en permanence dans lesdits signaux a debruiter (u(t)), en recherchant /Vtrames successives, 
de p echantillons chacune, ayant les caracteristiques attendues d'un bruit, en stockant les NxP echantillons cor- 
respondents pour constituer ledit modele de bruit, et en reiterant la recherche pour trouver un nouveau modele 
de bruit et stocker le nouveau modele en remplacement du precedent ou conserver le modele precedent selon 
les caracteristiques respectives des deux modeles. 

Application du procede selon Tune quelconque des revendications precedentes au debruitage de signaux de parole 
bruites (u(t)). 

Application du procede selon la revendication 8, caracterise en ce que la duree desdites frames est comprise dans 
la gamme 10 a 20 ms. 
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